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Abstract 

We consider tilings and packings of H d by integral translates of cubes [0,2[ d , which are 
4Z d -periodic. Such cube packings can be described by cliques of an associated graph, which 
allow us to classify them in dimension d < 4. For higher dimension, we use random methods 
for generating some examples. 

Such a cube packing is called non- extendible if we cannot insert a cube in the complement 
of the packing. In dimension 3, there is a unique non-extendible cube packing with 4 cubes. 
We prove that d-dimensional cube packings with more than 1 d — 3 cubes can be extended to 
cube tilings. We also give a lower bound on the number N of cubes of non-extendible cube 
packings. 

Given such a cube packing and z G Z rf , we denote by N z the number of cubes inside the 
4-cube z + [0, 4[ d and call second moment the average of N%. We prove that the regular tiling 
by cubes has maximal second moment and give a lower bound on the second moment of a 
cube packing in terms of its density and dimension. 

1 Introduction 

A general cube tiling is a tiling of by translates of the hypercube [0, 2[ d , which we call a 2-cube. 
A special cube tiling is a tiling of R d by integral translates of the hypercube [0,2[ d , which are 
4Z d -periodic. An example of such a tiling is the regular cube tiling of R d by cubes of the form 
z+ [0,2[ d with z G 2Z d . 

In dimension 1, there is only one type of special cube tiling, while in dimension 2, two following 
types of special cube tilings exist: 



The Keller's cube tiling conjecture (see [Kc30 ) asserts that any tiling of R d by translates of a unit 
cube admits at least one face-to-face adjacency. It is proved in |Sza86j that if this conjecture has a 
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Figure 1: The unique non-extendible cube packing in dimension 3 



counter example, then there is another counter example, which is also a special cube tiling. Using 
this, the Keller conjecture was solved negatively for d > 10 in [LaShQ2] and d > 8 in |McKa02j 
(note that the conjecture is proved to be true for d < 6 in |Pe40j ). Hence, special cube tilings, while 
seemingly limited objects have a lot of combinatorial possibilities. In the rest of this paper cube 
tiling stands for special cube tilings and iV is the number of orbits of cubes under the translation 
group 4Z d . Another equivalent viewpoint is to say that we are doing tilings of the torus R d /4Z d 
and N is then the number of cubes in this torus. 

A cube packing is a 4Z d -periodic set of integral translates of the 2-cube, such that any two 
cubes are non-intersecting. In dimension d > 3, there exist cube packings, called non-extendible, 
which cannot be extended to a tiling of the space (this first appear in |La00j ). In dimension 3 this 
non-extendible packing is unique (see Figure and it is the source of much of the inspiration of 
this paper. 

In Section |2l following [LaSh92j, we present a translation of the packing and tiling problems 
into clique problems in graphs. Explicit methods, in GAP, are used up to d = 4. For d > 5, we 
use various random methods, in Fortran 90 and C++ for generating random cube packings. 

Denote by f(d) the smallest number of cubes, which form a non-extendible cube packing. In 
Section El we give some lower and upper bounds on the value of f(d). In [DIP05 it is proved that 
any cube packing of [0, 4[ d by cubes [0, 2[ d is extendible to a 4Z d -periodic cube tiling of R d . 

If CV is a cube packing, denote by hole{CV) and call hole, its complement H d — CV. We prove 
that if a cube packing has more than 2 d — 3 cubes, then it is extendible to a tiling, i.e. that holes 
of volume at most 3 are tillable. We also obtain some conjecture on nonfillable holes of volume at 
most 7. 

Given a cube packing CV, the counting function N Z (CV) is defined as the number of cubes of 
CV contained in 2+ [0, 4[ d . We study its second moment in Section 0] We prove, that the highest 
second moment for tilings is attained for the regular cube tiling and give a lower bound for the 
second moment of cube packings, in terms of its dimension d and number of cubes N. 

2 Algorithm for generating cube packings 

Every 2-cube of a <i-dimensional cube packing is equivalent under 4Z d to a cube with center in 
{0, 1, 2, 3} d . Two 2-cubes of centers x and x' do not overlap if and only if there exist a coordinate 
i, such that \xi — x[\ =2. So, one consider the graphs Gd (introduced in |LaSh92]) with vertex-set 
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{0, l,2,3} d and two vertices being adjacent if and only if their associated cubes do not overlap. 
Cube packings correspond to cliques of Gd', they are non-extendible if and only if the cliques are 
not included in larger cliques. Cube tilings correspond to cliques of size 2 d . 

For a given d, the graph Gd has a finite number of vertices and a symmetry group Sym(Gd) 
of size d\.8 d . Hence, it is theoretically possible to do the enumeration of the cliques of Gd- The 
algorithm consists of using the set of all cliques with N vertices, considering all possibilities of 
extension, and then reducing by isomorphism using SymiGd) (the actual computation was done in 
GAP, see |GAPj ). The group SymiGd) is presented as a permutation group in GAP and the cliques 
as subsets of {1, ... , 4 d }. GAP uses backtrack search for testing if two subsets are equivalent under 
Sym(Gd) and is hence, very efficient even for large values of d. This enumeration is, in practice, 
possible only for d < 4 due to the huge number of cliques that appear. 

For d = 2, one finds only two non-extendible cliques of 4 vertices, i.e. two cube tilings. For 
d = 3, there is a unique non-extendible clique with 4 vertices, while there are 9 orbits of non- 
extendible cliques with 8 vertices (i.e. cube tilings). For d = 4, the computations are still possible 
and one finds the following results with N being the number of vertices of the clique. 



N 


1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 


#nonext. orbit cliques 


38 6 24 71 744 



Suppose that we have a cube tiling with two cube centers x and x' , satisfying to x' = x + 2ej with 
ej being the i-th unit vector, i.e. they have a face-to-face adjacency. If one replace x, x' by x + e iy 
x' + Ci and leave other centers unchanged, then one obtains another cube tiling, which we call the 
flip of the original cube tiling. The enumeration strategy is then the following: take as initial list 
of orbits the orbit of the regular cube tiling. For every orbit of cube packing, compute all possible 
pairs {x, x'}, which allow to create a new cube tiling. If the corresponding orbits of cube tilings 
are new, then we insert them into the list of orbits. Given a dimension d, consider the graph Cod, 
whose vertex-set consists of all orbits of cube tilings and put an edge between two orbits if one is 
obtained from the other by a flipping. The above algorithm consists of computing the connected 
component of the regular cube tiling in Cod- Since the Keller conjecture is false in dimension 
d > 8, we know that in those dimensions there are some isolated vertices in the graph and so, the 
above algorithm does not work. However, the graph Cod is connected for d < 4, i.e. any two cube 
tilings in those dimensions can be obtained by a sequence of flipping. It is an interesting question 
to decide, in which dimension d the graph Cod is connected; the only remaining unsolved cases are 
d = 5,6, 7. 

For dimensions d > 5, two above enumerative methods cannot work since there are too much 
possibilities. Hence, we used random methods. The random packing consist of selecting points, at 
random, on {0, 1, 2, 3} d , so that the corresponding 2-cubes do not intersect, until one cannot do 
this any more. This random packing algorithm creates non-extendible cube packings. 

The actual algorithm for creating non-extendible cube packings is as follows: the list L of 
selected cubes is, initially, empty. One select at random elements of {0, 1, 2, 3} d and keep them if 
they are adjacent to preceding elements of L. Of course not every trial works and as the space 
becomes more and more filled, the number of random generation needed to get a non-overlapping 
cube increase. When this number has reached a certain level, we go to a second stage: enumerate 
all possible extensions of the clique, that we have, and work in this list by eliminating elements 
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of it after choices are made. This algorithm has the advantage of enumerating the set {0, 1, 2, 3} d 
only one time and is hence, relatively fast. 

If one wants to find some packings with low density, then the above strategy is not necessarily 
the best. The greedy algorithm consists of keeping all 4 d elements in memory and at every step 
generate, say 20 elements and keep the one which cover the largest part of the remaining space. 

Another possibility is what we call Metropolis algorithm (see |Liu01| ): we take an non-extendible 
cube packing, remove a few cubes and rerun a random generation from the remaining cubes. If 
obtained packing is better than the preceding one, or not worse than a specified upper bound, then 
we keep it; otherwise, we rerun the algorithm. This strategy allows to make a random walk in the 
space of non-extendible cube packings and is based on the assumption, that the best non-extendible 
cube packings are not far from other, less good, non-extendible cube packings. 

3 Non-Extendible cube packings 

In dimension 1 or 2, any cube packing is extendible to a cube tiling. The exhaustive enumeration 
methods of the preceding section show that in dimension 3, there is a unique non-extendible cube 
packing. The set of its centers is, up to a symmetry of G 3 : 

{(0,0,0), (3, 2, 3), (2, 1,1), (1,3, 2)}. 

and its corresponding drawing is shown on Figure Its space group symmetry is P4(l)32, which 
is a chiral group. 

We first concentrate on the problem of finding non-extendible cube packings with the smallest 
number f(d) of cubes. From Sectional we know that /(l) = 2, f{2) = 4, /(3) = 4 and /(4) = 8. 

Lemma 1 For any n,m > 1, one has the inequality f(n + m) < f(n)f(m). 

Proof. Let Pa and Pb be non-extendible cube packings of R n and R m with f(n) and f(m) cubes, 
respectively. Let a k = (a k , a\, . . . , a k ) and b l = (b[, b\, . . . , b l m ) with 1 < k < f(n), 1 < I < f(m) 
be the centers of the 2-cubes from Pa and Pb- 

Define P to be the set of 2-cubes C kl with centers c kl = (a k , a^, ■ ■ ■ , a k ,b[,b l 2 , ■ ■ ■ ,b l m ) for 
1 < k < f{n) and 1 < I < f(m). The size of P is f{n)f(m) and it is easy to check that P is a 
packing. 

Take a cube D with center d = (di, e^, • • • , d n+m ). The vector (di, . . . , d n ) overlaps with a 
2-cube, say A k ° in Pa, while the vector (d n+ i, . . . ,d n+m ) overlaps with a 2-cube, say B l ° in Pb- 
Clearly, D overlaps with C k ° l ° and P is non-extendible. □ 

Since, /(3) = 4, one has /(6) < 16. 

A blocking set is a set {v>} of vectors in {0, 1, 2, 3} d , such that for every other vector v, there 
exist a j such that the 2-cubes of center v> and v overlap. A priori, the 2-cubes corresponding to 
the vector set {v>} can overlap; so, one has obviously h(d) < f(d). 

It is easy to see that h{2) = 3 and that any blocking sets of size 3 belong to one of two following 
orbits: 
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A slightly more complicated computation shows that h(3) = 4 and that any blocking set of size 4 
belong to one of three following orbits: 

{(0,0,0), (1,1,1), (2, 2, 2), (3, 3, 3)}, 
{(0,0,0), (1,1,1), (2, 2, 3), (3, 3, 2)}, 
{(0,0,0), (3, 2, 3), (2, 1,1), (1,3, 2)} 

Lemma 2 Let N satisfy the inequality < h(d), then one has h(d+ 1) > N. 

Proof. First h(d) > N if and only if, for any set P of N 2-cubes, there exists a 2-cube D, which 
does not overlap with any 2-cube from P. 

Let P be a set of N 2-cubes in torus T d+1 . Then at least 1"^] centers of them have Xd+i = t, 
for some t e {0,1,2,3}. Let us define another set P' of vectors by removing those vectors and 
the d + 1-th coordinate for the remaining vectors. Then P' consists of at most iV — |~~] = IW"] 
2-cubes. But |_^J < h(d); so, there exists a 2-cube C with center c = (c%,C2, ■ ■ ■ , q), which do 
not overlap with any 2-cube in P'. But then the 2-cube with center (c\, C2, . . . , c^, t + 2) does not 
overlap with any 2-cube from P. □ 



Theorem 1 For any d > 1, one has h(d + 1) > 



Ah(d)-1 
3 



+ 1. 



Proof. Let N 



4h(d)-l 



then it holds: 



3N 



4h(d)-l 



< 



Ah(d) - 1 



< h(d). 



And, from Lemma El we have, that h(d + 1) > N. □ 
Theorem Q does not allow to find an asymptotically better lower bound on f(d) than the trivial 
lower bound |~(§) d ]. Note that using Lemma [T] one proves easily that the limit (3 = lim^oo ln ^ 
exists. This limit satisfies to | < e^ < \fi. The upper bound following from Lemma ^ and 
/(3) = 4. The determination of (3 is open. 

Proposition 1 One has h(A) = 7. 

Proof. The following set of center coordinate proves that h(4) < 7. 

{(0, 0, 0, 0), (1, 1, 1, 1), (2, 2, 2, 2), (3, 3, 3, 3), (0, 0, 1, 1), (1, 1, 2, 2), (2, 2, 3, 3)} 

From Theorem Q we have h(4) > 6. Assume that h(4) = 6 and take a blocking set of six 2-cubes 
with centers a* = (a\, a\, a l 3 , a\), 1 < i < 6. 
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If three vectors a 11 , a 12 , a 23 have equal coordinate j, then by a reasoning similar to Lemma |2l 
one finds a vector which does not overlap with those vectors. 

So, the above situation does not occur and for every coordinate j, there exist two pairs {a ?1 , a 12 }, 
{a* 3 , a 14 }, which have equal j coordinates. 

We have two pairs A and B in first column. Take a pair A' in second column and assume that 
it does not intersect with A. Denote by P' the set of vectors obtained by removing the vector 
corresponding to the sets A and A' and the first and second coordinate of the remaining vectors. 
P' is a set of two vectors in dimension 2; hence, it is not blocking. So, we can find a 2-cube, which 
does not overlap with P. So, any of six pairs from three other columns must intersect with A and 
B. 

But we have only 4 different ways to intersect A and B. So, two pairs from column 2 — 4 are 
equal. But, if two pairs are equal, then they do not intersect, which is impossible. So, h(4) > 6. □ 
Theorem ^ and Proposition imply the following inequalities: 

/(5) > h(B) > 10 and /(6) > h(6) > 14 

By running extensive random computation we found more than 140000 non-extendible cube 
packings in dimension 5 with 12 cubes; they belong to 203 orbits. Hence, it seems reasonable 
for us to conjecture that in fact /(5) = 12 and that the number of orbits of non-extendible cube 
packings with 12 cubes is "small", i.e. a few hundreds. 

But dimension 6 is already very different. We know that /(6) < 16 but we are unable to find 
by random methods a single non-extendible cube packing with less than 20 cubes. 

We now consider cube packing with high density. 

Take a cube packing of R d with center set {x fc }, 1 < k < N. Select a coordinate % and an 
index j and form a cube packing of R^" 1 , called induced cube packing on layer j, by selecting all 
x fc with x\ = j,j + 1 (mod 4) and then creating the vector (a;*, . . . , . . . , x%). 

Lemma 3 If CV is a cube packing with 2 d — S cubes, then its induced cube packings have at least 
2 d - x - S cubes. 

Proof. Select a coordinate and denote by rij the number of 2-cubes of CV, with x% = j. One has 
no + n\ + n 2 + n 3 = 2 d - 5. 

The number of 2-cubes of the induced cube packing on layer j is yj = rij + nj + i. One writes 
yj = 2 d ~ l - Sj with 5j > 0, since the induced cube packing is a packing. Clearly, one has 
o~o + Sx + 5 2 + 5 3 = 26. 

We have rij +nj +1 = 2 d ~ 1 — 5f, so, one gets, by subtracting rij — rij +2 = Sj + i — Sj, which implies: 

5 — 5i + 5 2 — S 3 = . 

Every vector A = (So, 5±, 5 2 , 5 3 ) G Z^_, satisfying the above relation, can be expressed in the form 
c (l, 0, 0, 1) + ci(l, 1, 0, 0) + c 2 (0, 1, 1, 0) + c 3 (0, 0, 1, 1) with Cj e Z + . This implies Sj = Cj + c j+ i < 

V'-, ^ ' □ 

Theorem 2 In dimension d, every cube packing with 2 d — S cubes for S = 1, 2, 3 can be extended 
to a cube tiling. 
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Proof. The proof is by induction on d. Take d > 4 and a cube packing CV with 2 d — 5 cubes and 
denote by hole(CV) its hole in R d . Let us consider the layering along the coordinate i. By Lemma 
the induced cube packings have 2 d ~ 1 — 8j cubes with 5j < 3. So, one can complete them to form 
a cube packing of R^ 1 . Denote by CC, = [0, 2[ i ~ 1 x [0, l[x [0, 2[ d-< the half of a 2-cube cut along the 
coordinate i. The induced cube packings are extendible by the induction hypothesis. This means 
that hole(CV) is the union of 25 cut cubes CCi. Denote by Aj = (do, 8%, 82, 83) the corresponding 
vector; by the analysis of Lemma |3] Aj = c (l, 0, 0, 1) + ci(l, 1, 0, 0) + c 2 (0, 1, 1, 0) + c 3 (0, 0, 1, 1) for 
some Ci G Z + with = 5. 

In the case 5 = 1, it is clear that the only set, which for any i can be written as v 1 ' 1 + CCi U 
v 2,i + CCi for some vectors v 1 '', v 2 '* is the 2-cube itself. 

If 5 = 2, then, clearly, the vector A; takes, up to isomorphism, one of three different forms: 
(1,1,1,1), (1,2, 1,0) or (2, 2, 0,0). 

Suppose that for a given i, the vector Aj contains the pattern (0, 1). This means that on one 
layer we have exactly one translate, say v + CCi, of CCi. Select any other coordinate i' , v + CCi is 
splitted in two parts, say v 1 + CC^ and v 2 + CC^ by the layers along the coordinate i'. Since, an 
adjacent layer is completely filled, this means that v 2 = v 1 ± ey. Hence, they form a cube and the 
cube packing is extendible. 

Suppose that for a given coordinate i, Aj = (x, x, 0, 0) with x = 2 or 3. The 0-th layer is 
filled with x translates of set CCi. Take another coordinate, say i' , and consider the partition of 
hole(CV) into translates of CCi>. By intersecting with the 0-th layer, one obtains 2x intersections. 
But since the third layer is full, it is necessary for the translate of CC^ to overlap only on 1-th 
layer. This means that they make a cube tiling. 

The above considerations settle cases (2,2,0,0) and (1,2,1,0). Now assume that for a given 
coordinate i, one has Aj = (1,1, 1,1). Assume also that the cube packing is non-extendible. Take 
one translate v + CCi on layer j in hole(CV). It is splitted in two parts by the translates of CCi'. 
Since we assume that the cube packing is non-extendible, one of these translates overlaps on layer 
j — 1 and the other one on layer j + 1. One obtains a unique stair structure as illustrated below 
in a two-dimensional section: 

































□ 













Now select another coordinate i" (since d > 4) and see that hole(CV) cannot be decomposed into 
translates of CCi". So, if 5 = 2, then all cube packings are extendible. 

If 5 = 3, then for a given coordinate i, one has clearly, up to isomorphism, Aj=(3, 3, 0, 0), 
(2, 1, 1, 2) or (2, 3, 1, 0). The cases (3, 3, 0, 0) and (2, 3, 1, 0) are extendible by the above analysis. 
Let us consider the case (2, 1, 1, 2) and assume that the cube packing is non-extendible. The 1-th 
and 2-th layers consist of translates of CCi, which we write as v 1 +CCj and v 2 + CCj. The translate 
v 2 + CCi is splitted in two by the translate of CC^ appearing in the decomposition of hole(CV) 
along coordinate i'. If those translates spilled only on the 0-th layer or 2-th layer, then one has a 
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cube, which is excluded. So, they spill on 0-th and 2-th layers. This implies that v 2 = v 1 ± e^. 
But this is impossible, since %' is arbitrary. So, the cube packing is extendible. □ 
Given a d- dimensional non-extendible cube packing with 2 d — 5, its lifting is a d+ 1- dimensional 
non-extendible cube packing obtained by adding a layer of cube tiling; the iteration of lifting is 
also called lifting. 

Conjecture 1 Take CV a non-extendible cube packing with 2 d — 5 cubes. On its hole we conjecture: 

1. If 5 = 4 then hole{CV) is obtained as the hole of the lifting of the unique non-extendible cube 
packing in dimension 3. 

2. The case 5 = 5 does not occur. 

3. If 5 = 6, then hole(CV) is obtained as the hole of the lifting of one of two non-extendible 
cube packing in dimension 4. 

4- If 5 = 7, then hole{CV) is obtained as the hole of the lifting of a non-extendible cube packing 
in dimension 4. 

This conjecture is supported by extensive numerical computations. We can obtain an infinity 
of non-extendible cube packings with 2 d — 8 cubes by doing layering of two (d — 1) -dimensional 
non-extendible cube packings with 2 d ~ 1 — 4 cubes. This phenomenon does not appear for non- 
extendible cube packings with 2 d — 9 cubes, but we are not able to state a reasonable conjecture 
for this case. 

4 The second moment 

Given a cube packing CV and z G Z d , N Z (CV) is defined as the number of 2-cubes of CV contained 
in z+ [0,A[ d - 

Given a 4Z d -periodic function /, its average is 

w) = ^ E /(*)■ 

^e{0,l,2,3} d 

We denote rrii(CV) the i-moment of CV, i.e. the average of N l z (CV). 
Theorem 3 Let CV be a cube packing with N cubes. One has: 

od 

mi (CV) = (—)N and m x {CV) + N(N - l)2~ d + 2- d d{2q(q - 1) + rq} < m 2 (CV) 
with N = 4q + r,0<r<3. 

Proof. Take N 2-cubes A 1 , . . . , A N with centers a 1 , . . . , a^. The 4-cube a+ [0,4[ d with corner 
(a±, . . . , aa) contains the 2-cube with center b = (&!,..., 6 d ) if and only if cij ^ hi for every i. Take 
all 4-cubes C 1: . . . , C 4 d. 
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Every 2-cube A 1 is contained in 3 d 4-cubes C^. Denote by rij the number of 2-cubes A 1 , 
contained in the 4-cube Cj. By definition, the first moment has the expression: 

k 

The second moment is equal to m 2 (CP) = -h Z~2k n l- Let *y ^ e ^ ne numbers of 4-cubes containing 
the 2-cubes A 1 and AK One has the relation: 

y-y t _ n k (n k - 1) 

l<i<j<N k=l 

which implies 4 d mi(CV) + 2 = 4 d m2(C'P). Let us denote by the number of equal coordi- 
nates of the centers a 1 and d> . Then one has 

tij = (l ) ^2 d >2 d + 2 d - 1 ^ 

The above inequality becomes an equality for = or 1. Summing over i and j one obtains 

t lJ >N(N-l)2 d - 1 + 2 d - 1 to 

l<i<j<N l<i<j<N 

Let us denote by Ri the number of equal pairs in column /. By definition, one has clearly: 

d 

Y to = Y Rl ■ 

l<i<j<N 1=1 

Let us fix a coordinate I and denote by d u the number of entries equal to u in column I. One has, 
obviously: 

Ri = Y^ d u {d u -\) ^ and ^ 4 = iV 

u=0 u=0 

The Euclidean division N = 4q + r and elementary optimization, with respect to the constraints, 
allow us to write: 

Rl > 2q(q -l)+rq 

The proof follows by combining all above elements. □ 
Note that the value of mi (CV) was already obtained in jDIP05j . For a fix d and N, we do not 

know which cube packing minimize the second moment. However, we can characterize which cube 

tilings have the highest second moment in Theorem 
Consider the following space of functions: 

/:{0,1,2,3}^R. 
Vx G {0, 1, 2, 3} d one has E x+{ o,iy f(?) = 1 and f(x) > 

It is easy to see that cube tilings correspond to (0, 1) vector in Q. Therefore, the problem of 
minimizing the second moment over cube tilings is an integer programming problem for a convex 
functional. 
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Theorem 4 The regular cubic tiling is the cube tiling with highest second moment. 
Proof. Given a function / £ Q, let us define 

M(f\(<r\ - I + f( x + e i) if ^ = or 2 
Mi{f){x) - | if = 1 or 3 . 

The function Mj(/) belongs to Q. Geometrically M»(/) is the cube packing obtained by merging 
two induced cube packing on coordinate i and layer and 2. We will prove E(N z (Mi(f)) 2 ) > 
E(N z (f) 2 ). Without loss of generality, one can assume, i — 1. 
The key inequality, used in computation below, is: 

(x + Xi + x 2 ) 2 + (xi + x 2 + x 3 ) 2 + (x 2 + x 3 + x ) 2 + (x 3 + X + Xi) 2 
< 2(x + xi + x 2 + x 3 ) 2 + (x Q + xi) 2 + (x 2 + x 3 ) 2 if Xi > 0. 

Define f Z2 {zi) = E« 2 e{Q,i,2}rf-i f( z i> z 2 + u 2) and obtain: 

4^(^(M X (/)) 2 ) = E f{ o,i,2,3} d (E ue{ o,i, 2 }^i(/)(^ + ^)) 2 

= E 2l =o E^efo.i^^- 1 (E« 1= o E U2 e{o,i,2}<i-i M 1 (f)(z 1 +u x ,z 2 + u 2 )) 2 

= E aa6 {0,l,2,3}"-i E zl =o(E' 1 =0 M l(^2)(2 ; l + «l)) 2 

= E^ow^^dSi^/^K)) 2 + (E« 1=0 ^(^i)) 2 + (E' 1=2 /, 2 (^i)) 2 } 

> E^eWiAS^ E' 1= o(E 2 1= o + ^i)) 2 = 4^(iV 2 (/) 2 ) 

Hence, using the operation Mi . . . M^, we can only increase the second moment. So, one gets: 

EiN^M, . . . M d {f)) 2 ) > E(N z (f) 2 ) for all / G Q . 

It is easy to see that M\ . . . Md(f) is the function with f(x) — 1 if x is a (0,2) vector and 0, 
otherwise; hence, it corresponds to a regular cube tiling. □ 
Note that it is easy to see that m 2 = for the regular cube tiling. 
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